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Abstract 

The Vlasov equation for the colhsionless evolution of the single-particle probability distribu- 
tion function (PDF) is a well-known example of coadjoint motion. Remarkably, the property 
of coadjoint motion survives the process of taking moments. That is, the evolution of the mo- 
ments of the Vlasov PDF is also coadjoint motion. We find that geodesic coadjoint motion of 
the Vlasov moments with respect to powers of the single- particle momentum admits singular 
(weak) solutions concentrated on embedded subspaces of physical space. The motion and 
interactions of these embedded subspaces are governed by canonical Hamiltonian equations 
for their geodesic evolution. 
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1 Introduction 

The Vlasov equation. The evolution of N identical particles in phase space with coordinates 
{qi,Pi) i = 1,2, ...,A^, may be described by an evolution equation for their joint probability 
distribution function. Integrating over all but one of the particle phase-space coordinates yields 
an evolution equation for the single-particle probability distribution function (PDF). This is the 
Vlasov equation. 

The solutions of the Vlasov equation reflect its heritage in particle dynamics, which may 
be reclaimed by writing its many-particle PDF as a product of delta functions in phase space. 
Any number of these delta functions may be integrated out until all that remains is the dy- 
namics of a single particle in the collective field of the others. In plasma physics, this collective 
field generates the total electromagnetic properties and the self-consistent equations obeyed by 
the single particle PDF are the Vlasov-Maxwell equations. In the electrostatic approximation, 
these become the Vlasov-Poisson equations, which govern the statistical distributions of parti- 
cle systems ranging from integrated circuits (MOSFETS, metal-oxide semiconductor field-effect 
transistors), to charged-particle beams, to the distribution of galaxies in the Universe. 

A class of singular solutions of the VP equations called the "cold plasma" solutions have a 
particularly beautiful experimental realization in the Malmb erg-Penning trap. In this experi- 
ment, the time average of the vertical motion closely parallels the Euler fluid equations. In fact, 
the cold plasma singular Vlasov-Poisson solution turns out to obey the equations of point-vortex 
dynamics in an incompressible ideal flow. This coincidence allows the discrete arrays of "vortex 
crystals" envisioned by J. J. Thomson for fluid vortices to be realized experimentally as solutions 
of the Vlasov-Poisson equations. For a survey of these experimental cold-plasma results see jS]. 



Vlasov moments. The Euler fluid equations arise by imposing a closure relation on the first 
three momentum moments, or p— moments of the Vlasov PDF /(p, t). The zero-th p— moment 
is the spatial density of particles. The first p— moment is the mean momentum and its ratio 
with the zero-th moment is the Eulerian fluid velocity. Introducing an expression for the fluid 
pressure in terms of the density and velocity closes the system of p— moment equations, which 
otherwise would possess a countably infinite number of dependent variables. 

The operation of taking p— moments preserves the geometric nature of Vlasov's equation. 
It's closure after the first p— moment results in Euler 's useful and beautiful theory of ideal 
fluids. As its primary geometric characteristic, Euler's fluid theory represents fluid flow as 
Hamiltonian geodesic motion on the space of smooth invertible maps acting on the flow domain 
and possessing smooth inverses. These smooth maps (called diffeomorphisms) act on the fluid 
reference configuration so as to move the fluid particles around in their container. And their 
smooth inverses recall the initial reference configuration (or label) for the fluid particle currently 
occupying any given position in space. Thus, the motion of all the fluid particles in a container 
is represented as a time-dependent curve in the infinite-dimensional group of diffeomorphisms. 
Moreover, this curve describing the sequential actions of the diffeomorphisms on the fluid domain 
is a special optimal curve that distills the fluid motion into a single statement. Namely, "A fluid 
moves to get out of its own way as efficiently as possible." Put more mathematically, fluid flow 
occurs along a curve in the diffeomorphism group which is a geodesic with respect to the metric 
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on its tangent space supplied by its kinetic energy. 

Given the beauty and utility of the solution behavior for Euler's equation for the first 
p— moment, one is intrigued to know more about the dynamics of the other moments of Vlasov's 
equation. Of course, the dynamics of the the p— moments of the Vlasov-Poisson equation is one 
of the mainstream subjects of plasma physics and space physics. 

Summary. This paper formulates the problem of Vlasov p— moments governed by quadratic 
Hamiltonians. This dynamics is a certain type of geodesic motion on the symplectomorphisms, 
rather than the diffeomorphisms. The symplectomorphisms are smooth invertible maps acting 
on the phase space and possessing smooth inverses. We shall consider the singular solutions of 
the geodesic dynamics of the Vlasov p— moments. Remarkably, these equations turn out to be 
related to integrable systems governing shallow water wave theory. In fact, when the Vlasov 
p— moment equations for geodesic motion on the symplectomorphisms are closed at the level 
of the first p— moment, their singular solutions are found to recover the peaked soliton of the 
integrable Camassa-Holm equation for shallow water waves j2]. 

Thus, geodesic symplectic dynamics of the Vlasov p— moments is found to possess singular 
solutions whose closure at the fluid level recovers the peakon solutions of shallow water theory. 
Being solitons, the peakons superpose and undergo elastic collisions in fully nonlinear interac- 
tions. The singular solutions for Vlasov p— moments presented here also superpose and interact 
nonlinear ly as coherent structures. 

The plan of the paper follows: 

Section |21 defines the Vlasov p— moment equations and formulates them as Hamiltonian sys- 
tem using the Kupershmidt-Manin Lie-Poisson bracket. This formulation identifies the 
moment equations as coadjoint motion on a dual Lie algebra q* , in any number of 
spatial dimensions. 

Section 131 derives variational formulations of the p— moment dynamics in both their Lagrangian 
and Hamiltonian forms. 

Section|l]formulates the problem of geodesic motion on the symplectomorphisms in terms of the 
Vlasov p— moments and identifies the singular solutions of this problem, whose support 
is concentrated on delta functions in position space. In a special case, the truncation 
of geodesic symplectic motion to geodesic diffeomorphic motion for the first p— moment 
recovers the singular solutions of the Camassa-Holm equation. 

Section El discusses how the singular p— moment solutions for geodesic symplectic motion are 
related to the cold plasma solutions. By symmetry under exchange of canonical momentum 
p and position q, the Vlasov g— moments are also found to admit singular (weak) solutions. 

2 Vlasov moment dynamics 

The Vlasov equation may be expressed as 



Here the canonical Poisson bracket [ • , • ] is defined for smooth functions on phase space with 
coordinates {q,p) and f{q,p,t) is the evolving Vlasov single-particle distribution function. The 
variational derivative Sh/5f is the single particle Hamiltonian. 



dt L ' 5 / J dp dq 6f 



dld_6h 
dq dp 6f 



adi 



■5h/5f 
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A functional g[f] of the Vlasov distribution / evolves according to 



dg ffSg df . , ffSg r Sh 

Tt = JJ sjm "'"'•'=]] If V-Sf 



dqdp 



f 



In this calculation boundary terms are neglected upon integrating by parts and the notation 
(( • , • )) is introduced for the pairing in phase space. The quantity {g , h} defined in terms of 
this pairing is the Lie-Poisson Vlasov (LPV) bracket [IF. This Hamiltonian evolution equation 
may also be expressed as 

I = {5, M = ((/, ^ds,/s/^)) = - {{^d*sH/sff, 

which defines the Lie-algebraic operations ad and ad* in this case in terms of the pairing 
on phase space (( • , • )): s* x s i— > M. Thus, the notation ad^^^^j / in (|2.1|) expresses coadjoint 
action of 6h/5f € 5 on / E s*, where 5 is the Lie algebra of single particle Hamiltonian vector 
fields and s* is its dual under L? pairing in phase space. This is the sense in which the Vlasov 
equation represents coadjoint motion on the symplectomorphisms. 

2.1 Dynamics of Vlasov Moments 

The phase space moments of the Vlasov distribution function are defined by 



gram = JJ f {q,P) Q^P"" dq dp . 

The moments gfhm are often used in treating the collisionless dynamics of plasmas and 
particle beams j^l. This is usually done by considering low order truncations of the potentially 
infinite sum over phase space moments, 

oo oo 

9 ^ ^ 0, mm9 m,m > h ^ ^ bnn9n,n ) 

rh,m=0 n,n=0 

with constants flmm and with m, m, n, n = 0, 1, . . . . If /i is the Hamiltonian, the sum over 
moments g evolves under the Vlasov dynamics according to the Poisson bracket relation 

dg °° 

■^ = {5,^} = X] <^mmbnn{'mm-nn)gfh+n-l,m+n-l- 

m,m,n,n=0 

The symplectic invariants associated with Hamiltonian flows of the moments were discov- 
ered and classified in ,12j . 

2.2 Dynamics of Vlasov p— Moments 

The momentum moments, or "p— moments," of the Vlasov function are defined as 



Amiq,t) = J p"" f{q,p,t)dp, m = 0,l. 
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That is, the p— moments are g— dependent integrals over p of the product of powers p™, m = 
0, 1, . . . , times the Vlasov solution f{q,p,t). We shall consider functional of these p— moments 
defined by, 

5 = XI (^rriiq) P"^ f dqdp = j ara{q) Am{q) dq =: ^ l^Ara , a, 

m=0 m=0 m=0 

OO OO ^ OO 

^ = J2h^-(^')P''fdQdp = J2hniQ)MQ)dq=-.^{An,f3n 



n=0" 



n=0" 



n=0 



where ( • , • ) is the pairing on position space. 

The functions am and with m, n = 0, 1, . . . are assumed to be suitably smooth and 
integrable against the Vlasov p— moments. To assure these properties, one may relate the 
p— moments to the previous sums of Vlasov moments by choosing 



am{q) = ^afhmq'^, Pn{q) = ^hifinq' 



m=0 



n=0 



For these choices of oim{q) and Pniq), the sums of p— moments will recover the full set of Vlasov 
moments. Thus, as long as the moments of the distribution f{q,p) continue to exist 
under the Vlasov evolution, one may assume that the dual variables am{q) and Pn{q) are smooth 
functions whose Taylor series expands the p— moments in the g,p— moments. These functions 
are dual to the p— moments Am{q) with m = 0, 1, . . . under the pairing (• , •) in the spatial 
variable q. In what follows we will assume homogeneous boundary conditions. This means, for 
example, that we will ignore boundary terms arising from integrations by parts. 

2.3 Poisson bracket for Vlasov p— moments 

The Poisson bracket among the p— moments is obtained from the LPV bracket via the following 
explicit calculation, 



OO „„ 

{g,h} = - X jj f[o^m{q)p'\Pn{q)p 



dqdp 



m,n=0 

OO 



m,n=0 

OO 



X // mamPn' -nPnam' fp"''^'' ^dqdp 



m,n=0 

OO 

X \A-m+n-l , ad/3„Q;, 



dq 



m,n=0 

OO 



X / IT'PnA'm+n-l + {m + n)Am+n-lPn ' Oim dq 

n 



m,n=0 

OO 



m,n=0 



OLr 
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where we have integrated by parts and introduced the notation ad and ad* for adjoint and 
coadjoint action, respectively. Upon recahing the dual relations 



the LPV bracket in terms of the p— moments may be expressed as 

dh d . ... d 6h 



{g,h}{{A}) = - J n-^-^Am+n-i + {m + n)Am+n-i 



m,n=0 



dq 5Ar, 



dq 





5h 1 


UAm 


SAnl 



— • ^ ^ ^^m+n— 1 5 

m,n=0 

This is the Kupershmidt-Manin Lie-Poisson (KMLP) bracket which is defined for functions 
on the dual of the Lie algebra with bracket 

This Lie algebra bracket inherits the Jacobi identity from its definition in terms of the canonical 
Hamiltonian vector fields. Thus, we have shown the 

Theorem 2.1 (Gibbons [7j) 

The operation of taking p— moments of Vlasov solutions is a Poisson map. It takes the LPV 
bracket describing the evolution of f{q,p) into the KMLP bracket, describing the evolution of 
the p— moments An{x). A result related to this, for the Benney hierarchy 11], was also noted by 
Lebedev and Manin 114^ . 

The evolution of a particular p— moment Amiq,t) is obtained from the KMLP bracket by 
dAm r . , / 6h d d 6h 

-^ = {A^,h} = + 

71=0 

The KMLP bracket among the p— moments is given by 

{A^,A^} = 

oq oq 

expressed as a differential operator acting to the right. This operation is skew-symmetric under 
the L^ pairing and the general KMLP bracket can then be written as |7] 



m,n=0 

SO that 



dq 



A T — 

dt -Z^^^-'^^^M,' 

n=0 
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2.4 Multidimensional treatment 

We now show that the KMLP bracket and the equations of motion may be written in three 
dimensions in multi-index notation. By writing p^^+i = ^^^^ checking that: 

i+j+k=n 



it is easy to see that the multidimensional treatment can be performed in terms of the quantities 
where cr G N^. Let A^- be defined as 

(q,t) =: Jp''f{(l,p,t)dp 

and consider functionals of the form 

g = ^ a^{(i) fj (q, p, i) dqdp =: ^ {A^, a„) 

h = J2ff^P (q)p''/ (q, P, t) dqdp =: ^ {A^, (3^) 
The ordinary LPV bracket leads to: 



{9M = -Y. jj {ci)p'',(3p{ci)pP]dqdp = 

= -EE/// "'"'"^r"^"'''''"'''^^ 



n d . , , , df3p 



EE 

o-,p i 



a„dq 



where the sum is extended to all cr, p G and we have introduced the notation, 

l-=:(0,... ,1, -,0) 



/'^ element 
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so that (Ij)^ = 5ji. 

The LPV bracket in terms of the p-moments may then be written as 



dt 



pGN3 j 

where the Lie bracket is now expressed as 





5h ] 







3 



d 6h d 6g 



Moreover the evolution of a particular p-moment A^j is obtained by 
dA^ 



dt 



{A.,h} = 

-EE 



6h d d 5h 



5Ap dqj ' 



J 



j dqj 5Ap 



and the KMLP bracket among the multi-dimensional moments is given in by 

{A., M = -E {^^^^^<^^P'^j + p^^'^^p-^j I") • 

Inserting the previous operator in this multi-dimensional KMLP bracket leads to 

{9,h}i{A}) = Y^J^JA.,Ap} '^^ ■ 



a,p 

and the corresponding evolution equation becomes 

dA^ ^ 

p 



6Ar 



dt 



6h 



Thus, in multi-index notation, the form of the Hamiltonian evolution under the KMLP bracket 
is essentially unchanged in going to higher dimensions. 

2.5 Applications of the KMLP bracket 

The KMLP bracket was derived in the context of Benney long waves, whose Hamiltonian is 

H2 = ^{A2 + Al). 



This leads to the moment equations 

dAn ^ dAn+i 
dt dq 



dAo 

+ nA„_i— ^ = 
dq 



derived by Benney |T| as a description of long waves on a shallow perfect fluid, with a free surface 
at y = h{q,t). In his interpretation, the An were vertical moments of the horizontal component 
of the velocity p{q, y, t): 



Ar,. = I p{q,y,tYdy. 

ly=0 
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The corresponding system of evolution equations for p{q, y, t) and h{q, t) is related by hodograph 
transformation, y = fi<l,p',t) dp', to the Vlasov equation 

^+ ^_^^ = 
dt dq dq dp 

The most important fact about the Benney hierarchy is that it is completely integrable. This 
fact emerges from the following observation. Upon defining a function X{q,p,t) by the principal 
value integral, 

\{q,p,t)=p + P j-dp , 

J^oo P-P 

it is straightforward to verify 14^ that 

dX dX dAo dX _ ^ 
dt ^ dq dq dp ' 

so that / and A are advected along the same characteristics. 

In higher dimensions, particularly n = 3, we may take the direct sum of the KMLP bracket, 
together with with the Poisson bracket for an electromagnetic field (in the Coulomb gauge) 
where the electric field E and magnetic vector potential A are canonically conjugate; then the 
Hamiltonian 



H 



MV 



2m 



1p - eA| 



„ n n 

/(p,q)d"pd"q+ / + -Y^^iAi, - Aj,,)' 

i=l j=l 



d"q 



yields the Maxwell- Vlasov (MV) equations for systems of interacting charged particles. For a 
discussion of the MV equations from a geometric viewpoint in the same spirit as the present 
approach, see For discussions of the Lie-algebraic approach to the control and steering of 
charged particle beams, see [5]. 



3 Variational principles and Hamilton-Poincare formulation 

In this section we show how the p— moment dynamics can be derived from Hamilton's principle 
both in the Hamilton-Poincare and Euler-Poincare forms. These variational principles are de- 
fined , respectively, on the dual Lie algebra g* containing the moments, and on the Lie algebra g 
itself. For further details about these dual variational formulations, see ^ and Summation 
over repeated indices is intended in this section 



3.1 Hamilton-Poincare hierarchy 

We begin with the Hamilton-Poincare principle for the p— moments written as 

6 dt{{An,[in)-H{{A}))={) 

Jti 

(where /3n £ s)- We shall prove that this leads to the same dynamics as found in the context 
of the KMLP bracket. To this purpose, we must define the n— th p— moment in terms of the 
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Vlasov distribution function. We check that 

ft 



= S f'dt{{An,Pn)-H{{A})) = 

Jti 

^£'*(M(/.P»M-(*/.f)»^ 

Now recall that any g = 6G/5f belonging to the Lie algebra s of the symplectomorphisms 
(which also contains the distribution function itself) may be expressed as 

SG 6G 



of OAm 



by the chain rule. Consequently, one finds the pairing relationship, 

Sf,(p-Pn-^4)\ = (sAn,(Pn 



Next, recall from the general theory that variations on a Lie group induce variations on its Lie 
algebra of the form 

Sw = u + [g, u] 

where u,w E s and u vanishes at the endpoints. Writing u = p^ijm then leads to 

r dt ({/, 5 (p^Pn))) = r dt ((/, {il + b"/3n, «]))) = 
Jti Jti 

= - I dt {J^Am, - (^n+m-1, [/^n, '7ml)) = 

J ti 

Consequently, the Hamilton-Poincare principle may be written entirely in terms of the moments 

as ^ 

SS = J^' dt!^ (5 An, (^Pn - 1^ j ^ - ( (im + &d*p^Am+n-l) , r?m) | = 

This expression produces the inverse Legendre transform 

^" " 5 An 

(holding for hyperregular Hamiltonians). It also yields the equations of motion 

dAm _ 1* 4 

- -aap^Ani+n-l 

which are valid for arbitrary variations 5Ajn and variations 5/3^ of the form 

SPm ='nm + ^Pn'Hm-n+l 

where the variations r)m satisfy vanishing endpoint conditions, 

Vm\t=ti ~ Vm\t=tj ~ 

Thus, the Hamilton-Poincare variational principle recovers the hierarchy of the evolution equa- 
tions derived in the previous section using the KMLP bracket. 



Gibbons, Holm & Tronci 



Singular geodesic Bows of Vlasov momentum moments 



11 



3.2 Euler-Poincare hierarchy 

The corresponding Lagrangian formulation of the Hamilton's principle now yields 



S £ L ({/?}) dt = £ lJL^5(5^)dt = 

I 5L 



t 

t 



llm—n+l 

)\dt 



= -r((Hi:+'^**^)'''")* 

upon using the expression previously found for the variations SPrn and relabeling indices appro- 
priately. The Euler-Poincare equations may then be written as 



d 5L 5L 



with the same constraints on the variations as in the previous paragraph. Applying the Legendre 

transformation 

A 

■^m — r- 

darn 

yields the Euler-Poincare equations (for hyperregular Lagrangians) . This again leads to the 
same hierarchy of equations derived earlier using the KMLP bracket. 

To summarize, the calculations in this section have proven the following result. 

Theorem 3.1 With the above notation and hypotheses of hyperregularity the following state- 
ments are equivalent: 

1. The Euler-Poincare Variational Principle. The curves (3n{t) o-re critical points of the 
action ^ 

5 r L ({/?}) dt = Q 
Jti 

for variations of the form 

SPm = 'nm + ad/3„r/m-n+l 

in which rjm vanishes at the endpoints 

'nm\t=ti — 'nm\t=tj — 

and the variations SA^ are arbitrary. 

2. The Lie-Poisson Variational Principle. The curves are critical points of the 
action 



5 j\{An,l3n)-H{{A}))dt = 

for variations of the form 
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where r]m satisfies endpoint conditions 

'>lin\t=ti — 'nm\t=tj — 

and where the variations 6An are arbitrary. 
3. The Euler-Poincare equations hold: 

d* -0 



4- The Lie-Poisson equations hold: 



Am = -a.d*sH A, 



SH -ri-m+n—l 

SAn 



For further details on the proof of this theorem we address the reader to . An analogous 
result is also valid in the multidimensional case with slight modifications. 

4 Quadratic Hamiltonians 

4.1 Geodesic motion 

We shall consider the problem of geodesic motion on the space of p— moments. For this, we 
define the Hamiltonian as the norm on the moment given by the following metric and inner 
product, 



h = -\\A\ 



MQ)Gnsiq,q')As{q')dqdq' (4.1) 



?i,s=0 



The metric GnsiQ^Q') is chosen to be positive definite, so it defines a norm for {A} £ q*. The 
corresponding geodesic equation with respect to this norm is found as in the previous section to 
be, 

= {Am-, h} = -22, \ nl3n-Q-Am+n~i + {m + n)Am+n-i-Q-l3nj (4.2) 

n=0 

with dual variables Pn ^ 9 defined by 

Xh OO ^ OO 

Pn = Yr=Y. Gns{q,q')As{q')dq' = ^Gns*As. (4.3) 
Thus, evolution under ()4.2() may be rewritten as coadjoint motion on g* 

-^ = {Am^h}=:-Y, ad^„ (4.4) 

n=0 

This system comprises an infinite system of nonlinear, nonlocal, coupled evolutionary equations 
for the p— moments. In this system, evolution of the m^^ moment is governed by the potentially 
infinite sum of contributions of the velocities [3n associated with n*^ moment sweeping the 
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(m+n— 1)*^ moment by coadjoint action. Moreover, by equation (|4.3jl . each of the /?„ potentially 
depends nonlocally on all of the moments. 

Equations (|4.1|) and (|4.3j) may be written in three dimensions in multi-index notation, as 
follows: the Hamiltonian is given by 

h=^\\A\f = ^Y1 [[ ^f^ (q-*) G^.u (q,q') A, (q', t) dqdq' 

so the dual variable is written as 

Pp = — = Y, Gpu (q,q') A, (q',t) dqdq' = ^ Gp, * 



4.2 Singular geodesic solutions 

Remarkably, in any number of spatial dimensions, the geodesic equation 1)4. 2|) possesses exact 
solutions which are singular; that is, they are supported on delta functions in g— space. 

Theorem 4.1 (Singular solution Ansatz for geodesic flows of Vlasov p— moments) 
Equation i4-^ admits singular solutions of the form 

N 

Aa{q,t) = / P/(q,t,a,)(5(q-Q,(t,a,))da,- (4.5) 

in which the integrals over coordinates Oj are performed over N embedded subspaces of the 
q— space and the parameters {Qj,Pj) satisfy canonical Hamiltonian equations in which the Hamil- 
tonian is the norm h in \4-l\ ) evaluated on the singular solution Ansatz i4.5\ ). 

In one dimension, the coordinates aj are absent and the singular solutions in (|4.5|) reduce to 

N 

As{q,t) = ^P/(g,t)<5(g-Q,(t)). (4.6) 
i=i 

In order to show this is a solution in one dimension, one checks that these singular solutions 
satisfy a system of partial differential equations in Hamiltonian form, whose Hamiltonian couples 
all the moments 

oo N 

= ^Y. ^ ^j'(Qj^^'^^^^PkiQkit),t)GUQj(.t),Qk{t)) 

n,s=0 j,k=l 

One forms the pairing of the coadjoint equation 

Am = — ^ ^^Gns*As^'m+n-l 
n,s 

with a sequence of smooth functions {ifm (q)}, so that: 

{A m+n— 1 5 

n,s 
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One expands each term and denotes Pj{t) := Pj{Qj,t): 

dqiPm iq)-^^[Pp{q,t)6iq-Qj)] 



^ / dqifrniq) 
3 '' 



Hq-Qj) 



Qpra _ 



dt 



dP"^ ~ 

-^Vm {Qj)+PrQjVm'{Qj 



Similarly expanding 

{Am+n-i,^^dG„s*AsV>m) = Yl j dq P^ Pp+''-' 5 (q - Q , 

n(pm' (Qj) Gns {Qj,Qk) -m(pm {Qj] 



niprn 'Gns Qk) " rrvp, 



dGns (q, Qk 



^ pspm+n-l 



dq 

dGns {Qj, Qk) 

dQj 



leads to 



.. dQj 
dt 



dPp 
dt 



n,s k 



-m 



n,s k 



m+n-1 dGns {Qj: Qk) 

dQj 



so that we finally obtain equations for Qj and Pj in canonical form, 



dQj dH, 



N 



dt dPj 



dt 



OH 



N 



dQj 



Remark about higher dimensions The singular solutions H4.5() with the integrals over 
coordinates Oj exist in higher dimensions. The higher dimensional singular solutions satisfy 
a system of canonical Hamiltonian integral-partial differential equations, instead of ordinary 
differential equations. 

5 Discussion 

5.1 Remarks about EPSymp and connections with EPDiff 

Importantly, geodesic motion for the p— moments is equivalent to geodesic motion for the Euler- 
Poincare equations on the symplectomorphisms (EPSymp) given by the following Hamiltonian 

H[f] = \ II f{q,P,t)g{q,p,q',p')f{q',p',t)dqdpdq'dp' (5.1) 
The equivalence with EPSymp emerges when the function Q is written as 

g{q,q',p,p')=Y,P''Gn^'n{q,q')p'"'. 
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Thus, whenever the metric G for EPSymp has a Taylor series, its solutions may be expressed in 
terms of the geodesic motion for the p— moments. 

Moreover the distribution function corresponding to the singular solutions for the moments 
is a particular case of the cold-plasma approximation, given by 

j 

where in our case a summation is introduced and p is written as a Lagrangian particle-like 
density: 

p,{q,t)=S{q-Qj{t)) 

To check this is a solution for the geodesic motion of the generating function, one repeats exactly 
the same procedure as for the moments, in order to find the following Hamiltonian equations 

dt ~ dP^ ^Z^"^^'^^'^' dt - dQ, 

where Pj = Pj o Qj denotes the composition of the two functions Pj and Qj. This procedure 
recovers single particle motion for density pj defined on a delta function. 

As we shall show below, these singular solutions of EPSymp are also solutions of the Euler- 
Poincare equations on the diffeomorphisms (EPDiff), provided one truncates to consider only 
first order moments ^U]. With this truncation, the singular solutions in the case of single-particle 
dynamics reduce in one dimension to the pulson solutions for EPDiff |2j. 

5.2 Exchanging variables in EPSymp 

One can show that exchanging the variables q ^ p in the single particle PDF leads to another 
nontrivial singular solution of EPSymp, which is different from those found previously. To see 
this, let / be given by 

f{q, P,t) = Y, ^{Q - Q3 (p. 0) - Pj (t)) 
j 

At this stage nothing has changed with respect to the previous solution since the generating 
function is symmetric with respect to q and p. However, inserting this expression in the definition 
of the m— th moment yields 

j 

which is quite different from the solutions found previously. One again obtains a canonical 
Hamiltonian structure for Pj and Qj. 

This second expression is an alternative parametrisation of the cold-plasma reduction above and 
it may be useful in situations where the composition Qj o Pj is more convenient than Pj o Qj . 

5.3 Remarks about truncations 

The problem presented by the coadjoint motion equation ()4.4() for geodesic evolution of p— moments 
under EPDiff needs further simplification. One simplification would be to truncate the (doubly) 
infinite set of equations in ()4.4() to a finite set. These moment dynamics may be truncated at 
any stage. For example, two alternatives could be 
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• by modifying the Hamiltonian, so that some moments are decoupled (or, perhaps so that 
some moments slaved to others). For example, one might also arrange for the metric Gns 

— to only couple nearest neighbors, or 

— to be diagonal, or 

— to be a multiple of the identity 

• by modifying the Lie-algebra in the KMLP bracket to vanish for weights m + n — 1 greater 
than a chosen cut-off value. 

5.4 Examples of simplifying truncations and specializations. 

For example, if we truncate the sums to m, n = 0,1,2 only, then equation 1)4. 4|) produces the 
coupled system of partial differential equations, 

^ = -ad}^Ao - ad^^Ai , 
dAi 

-Qf = -adft^o - ad^^^i - ad;^^2 , 
dA 

= -ad^^Ai -ad^^yia. 
Expanding now the expression of the coadjoint operation 

adp^Ak+h-i = (k + h) A^+h-idgPh + hPhdqAk+h-i 



and relabeling 

one calculates 



dAp 

dt 
dAi 

"df 
dM 
dt 



ad^^^fc = (A; + 1) AkdgPh + hPhdqAk 

-d,{Aopi)-2Aidgp2-2/32dgAi 

-AodgPo - 2Aidq(3i - (5idgAi - 3^25g/?2 - 2/325^^2 

-2Aidq(3p - 3A2d,pi - /3idgA2 



We specialize to the case that each velocity depends only on its corresponding moment, so 
that Ps = G * Ag, s = 0,1,2. If we further specialize by setting and A2 initially to zero, then 
these three equations reduce to the single equation 

dA^ 

= -PldgAl- 2AldgPl. 

Finally, if we assume that G in the convolution Pi = G * ^1 is the Green's function for the 
operator relation 

Ai = {l- a^dl)pi 

for a constant lengthscale a, then the evolution equation for Ai reduces to the integrable 
Camassa-Holm (CH) equation [2 in the absence of linear dispersion. This is the one-dimensional 
EPDiff equation, which has singular (peakon) solutions. Thus, after these various specializa- 
tions of the EPDiff p— moment equations, one finds the integrable CH peakon equation as a 
specialization of the coadjoint moment dynamics of equation (|4.4|) . 



Gibbons, Holm & Tronci Singular geodesic Bows of Vlasov momentum moments 17 

That such a strong restriction of the p— moment system leads to such an interesting spe- 
cial case bodes well for future investigations of the EPSymp p— moment equations. Further 
specializations and truncations of these equations will be explored elsewhere. Before closing, 
we mention one or two other open questions about the solution behavior of the p— moments of 
EPSymp. 

5.5 Open questions for future work 

Several open questions remain for future work. The first of these is whether the singular solutions 
found here will emerge spontaneously in EPSymp dynamics from a smooth initial Vlasov PDF. 
This spontaneous emergence of the singular solutions does occur for EPDiff. Namely, one sees 
the singular solutions of EPDiff emerging from any confined initial distribution of the dual 
variable. (The dual variable is fluid velocity in the case of EPDiff). In fact, integrability 
of EPDiff in one dimension by the inverse scattering transform shows that only the singular 
solutions (peakons) are allowed to emerge from any confined initial distribution in that case 
P]. In higher dimensions, numerical simulations of EPDiff show that again only the singular 
solutions emerge from confined initial distributions. In contrast, the point vortex solutions 
of Euler's fluid equations (which are isomorphic to the cold plasma singular solutions of the 
Vlasov Poisson equation) while comprising an invariant manifold of singular solutions, do not 
spontaneously emerge from smooth initial conditions in Euler fluid dynamics. Nonetheless, 
something quite analogous to the singular solutions is seen experimentally for cold plasma in a 
Malmberg- Penning trap [S]. Therefore, one may ask which outcome will prevail for the singular 
solutions of EPSymp. Will they emerge from a confined smooth initial distribution, or will they 
only exist as an invariant manifold for special initial conditions? Of course, the interactions of 
these singular solutions for various metrics and the properties of their collective dynamics is a 
question for future work. 

Geometric questions also remain to be addressed. In geometric fluid dynamics, Arnold and 
Khesin formulate the problem of symplecto-hydrodynamics, the symplectic counterpart of 
ordinary ideal hydrodynamics on the special diffeomorphisms SDiff. In this regard, the work of 
Eliashberg and Ratiu j2()j showed that dynamics on the symplectic group radically differs from 
ordinary hydrodynamics, mainly because the diameter of Symp(M) is infinite, whenever M is 
a compact exact symplectic manifold with a boundary. Of course, the presence of boundaries 
is important in fluid dynamics. However, generalizing a result by Shnirelman [^, Arnold and 
Khesin point out that the diameter of SDiff(M) is finite for any compact simply connected 
Riemannian manifold M of dimension greater than two. 

In the case under discussion here, the situation again differs from that envisioned by Eliash- 
berg and Ratiu. The EPSymp Hamiltonian (|5.H) determines geodesic motion on Symp(T*M'^), 
which may be regarded as the restriction of the Diff(T*R^) group, so that the Liouville volume 
is preserved. The main difference in our case is that M = T*M^ is not compact, so one of the 
conditions for the Eliashberg-Ratiu result does not hold. Thus, one may ask, what are the 
geometric properties of Symp acting on a symplectic manifold which is not compact? What 
remarkable differences between Symp and SDiff remain to be found in such a situation? 

Yet another interesting case occurs when the particles undergoing Vlasov dynamics are con- 
fined in a certain region of position space. In this situation, again the phase space is not compact, 
since the momentum may be unlimited. The dynamics on a bounded spatial domain descends 
from that on the unbounded cotangent bundle upon taking the p-moments of the Hamiltonian 
vector field. Thus, in this topological sense p-moments and q-moments are not equivalent. In 
the present work, this distinction has been ignored by assuming either homogeneous or periodic 
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boundary conditions. 
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